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1. Introduction 

This paper purports to present an explicit construction of the representations of the 
1-D Schrodinger group, Schi, in the basis of the state vectors of the harmonic oscillator. 
In physical terms, this entails for instance, the determination of transition amplitudes 
for oscillators subjected to external forces and varying frequencies. The ensuing results 
are of relevance to many areas of physics. Interestingly, our study also contributes 
to the mathematical theory of matrix orthogonal polynomials that is currently being 
actively developped. Indeed, as will be seen, the group-valued matrix elements of 
the Schi-representation are expressed in terms of such matrix orthogonal polynomials 
(MOPs). Furthermore, the natural Lie-algebraic setting of our presentation allows to 
systematically obtain various structural formulas that characterize these MOPs. Much 
efforts are currently being deployed to obtain MOPs with such properties and their 
symmetry analysis has been rather limited. It is hence remarkable that our study 
provides in a most simple physical context, such a case-study of matrix orthogonal 
polynomials. Let us now give some background for our analysis. 



1.1. The Schrodinger group Schi and its algebra 

The Schrodinger group Schi is defined as the (maximal) group of symmetry 
transformations of the free one- dimensional Schrodinger equation: 

'dt~ 2 9x2- 

This group was obtained by Niederer in 1972, who found that the familiar Galilei 
transformations 

{t,x) ^ {t + b,x + vt + a) (L2) 
are supplemented by 

dilations : (t, x) — > (d^t, dx) (1.3a) 
and 

expansions : (t,x) — ^ ( , — ^ — ) . (1.36) 

\ 1 + at 1 + at J 

The infinitesimal symmetry generators T, P, G, K., and E respectively associated 
to time and space translations, Galilei boosts, dilations, expansions and phase 
transformations are given by 



T =idt, D = i(2tdt + xd^ + ]^ 



K = -i(t'd, + txd,,+ ^^-'^ (1.4) 

P = —idx, G = —itdx + X, E = 1. 

They are easily seen to form an algebra which is isomorphic to the semi-direct sum of 
sl(2,R) and the Heisenberg algebra : sl(2,]R) ix Hi. The generators P, G and E form a 
basis for the Hi subalgebra, while if, D, and K close onto sl(2,R) under commutation. 
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1.2. Schi and harmonic oscillators 

A year later, Niederer [2] again, studied the maximal symmetry group of the Schrodinger 
equation 



dt 



1^2^ 1 2, 



1.5) 



with a harmonic potential to find that this equation admits an invariance group which 
is isomorphic to Schi. 

Since this is the realization of Schi that we shall use in this paper, we shall record 
for reference, basic facts about quantum-mechanical harmonic oscillators. As is well 
known, they are described in terms of the anhilation and creation operators 



J_f d_ 

72 V dx 



V2 \ dx) 



that satisfy the Heisenberg algebra commutation relation 

[a, a"*"] = 1 
and in terms of which 



;i.6) 

1.7) 



H = - 



1 92 1 , 

1 T 

2dx^ 2 



a a + 



1 



We shall use throughout the standard representation that has for basis the normalized 
number vector states {|n), n = 0, 1, 2, . . .}, {m\n) = 6mn, in which 

a\n) = \pn |n — 1), a^\ri) = \/nW\n + 1), if |n) = ^ra + (1.9) 

In the position representation, the wave functions are expressed in terms of the Hermite 
polynomials Hn{x) (see Appendix A.l ): 

1 



^n{x) = {x\n) 



X 



1.10) 



We shall also make use of the coherent states \z), z G C defined as the eigenstates of 
the annihilation operator a. These have the following properties: 



1^) = e"" |0), 



a\z) = z\z), 



n=0 



\z) = \z + f3), 



{n\z) 



f{a)\z) = f{z)\z), 



;i.lla) 

(1.116) 
(1.11c) 



where / is an arbitrary (smooth) function in ( ll.ll&p . 

Let us now give an indication as to why the invariance group of fll.Sp is isomorphic 
to Schi. The symmetry transformations of space and time are obviously different in this 
case than those that map solutions of the free equation onto themselves, it is however 
easy to see that their generators span the same abstract algebra, that is sl(2,R) x Hi. 

We already noted that a, and / form a basis for the Heisenberg algebra Hi. 
Similarly the bilinears a^, a'^a and (a"*")^ realize the sl(2, M) commutation rules. Observe 
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that the sets of state {\2n),n = 0, 1, . . .} and {\2n + 1), n = 0, 1, . . .} with even and odd 
quanta, will separately transform onto themselves under the action of those bilinears; 
as a matter of fact, they span two different irreducible discrete representation series of 
sl(2,M). It is straightforward to check that together the linears and bilinears in a and 
realize the semi-direct sum sl(2,]R) x Hi. 
Consider the symmetry generators (11.41) of the free Schrodinger equation at t = 

P(0) = -id^ G{0) = X, E{0) = /, 

1 1 (1-12) 

m = --dl D{0) = xd^ + -, K{0) = . 

We know they form a sl(2, M) x Hi algebra. They also clearly span the set of linears 

and bilinears in a and a"*" at since 

i 1 

P(0) = -(a-a+), G'(O) = ^(a + a+), E(0) = aa+ - a+a, 

V 2 V 2 

T{0) = -^{a^ + {a+f-aa+ -a+a), D{0) = ^{a^ - {a+f), (1.13) 

K{0) = -^(a^ + {a+f + aa+ + a+a). 

If we take these operators generically denoted by X{0) and evolve them in time using 
the harmonic oscillator Hamiltonian H given in (11.81) (rather than the free Hamiltonian) 
according to 

X{t)= e'"'X{0)e-'"'; (1.14) 
this will obviously yield constants of motion, 

^(t) = ^(t)-z[i7,X(t)] = 0, (1.15) 

that is, symmetry generators. Since the time evolution (11.141) is given by a unitary 
automorphism, the time-dependent symmetry generators will satisfy at arbitrary t 
the same commutation relations than at t = 0, namely those of sl(2,M) k Hi. The 
symmetry transformation of (t, x) that leave (ll.lSp invariant can then be obtained 
by exponentiating these generators. This explains why the symmetry groups of the 
Schrodinger equation for the free particle and for the harmonic oscillator are structurally 
identical. 



1.3. Representations of the Schi on oscillator states: the central object 

The projective representations of the Schrodinger group in n dimensions Sch„ have 
been classified by Perroud [3]. More recently, an infinite-dimensional extension of the 
Schrodinger algebra known as the Schrodinger- Virasoro algebra has been constructed 
and studied [Sj HJ [5]. This extended algebra has been used in particular, to analyze 
anisotropic critical systems and to formulate a non-relativistic version of the AdS / CFT 
correspondance. The book of Roger and Unterberger [6] covers the mathematics of these 
developments; moreover, a short review of the many physical situations (old and new) 
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exhibiting Schrodinger symmetry, with references to the hterature, is provided in the 
introduction. 

In the following, we shall be interested in the explicit determination of the matrix 
elements for the unitary representation of Schi in the basis {|n), n = 0, 1, 2, . . .} of the 
oscillator quantum states. A priori, 6 parameters are needed to specify a general Schi 
group element. However, since the factor e*^^ ■ e*"^ is diagonal and hence trivial in 
the oscillator basis, it will not be included. With this understanding, we shall use the 
following parametrization for a general group element of Schi: 

S{v, W) = e^™-^'*'') ^{wa^-W{a+f )/2 ^ g^j^^ (^_^g^ 

where v and w are complex numbers (equivalent to four real parameters) and v and w 
their complex conjugate. Note that S' is a unitary operator : = S~^. 
Our central objects of study will hence be the matrix elements 

^Pn,k = {k\S\n) = {k\ e^-"^^ e^"'"'-^^"^)')/^). (1.17) 

As already mentioned, we shall show that these quantities are naturally expressed 
in terms of 2 x 2 matrix orthogonal polynomials Vn{k) of the discrete variable k. The 
reader will also have realized the relations with the theory of coherent and squeezed 
states [7] that our analysis will entail. 

1.4- Matrix orthogonal polynomials 

Introduced by Krein some 60 years ago, the theory of matrix orthogonal polynomials 
(MOPs) had not much progressed until a few years ago when Duran, Griinbaum and 
coworkers gave it a significant impetus with the discovery of families of MOPs that 
satisfy second-order differential equations. A survey of these advances can be found 
in [H] and our study adds to this corpus of results. Let us also signal Ref. |H], were 
the relation between polynomials obeying higher-order recurrence relations and matrix 
orthogonal polynomials was identified. 

We offer below basic definitions and properties of matrix orthogonal polynomials. 
(The reader is referred to [TO] for a broad exposition.) Let ('P„)„ be a sequence of 
polynomials with N x N matrix coefficients (or equivalently, a sequence of N x N 
matrices with polynomials of degree n or less as entries). In the case where these Vn 
depend on a continuous variable t, they will be orthogonal if there exists a suitably 
defined weight matrix W{t) such that 

J Vnit) m{t)Vi{1:) = hJ^^Ul)NxN. (1.18) 

Like in the scalar case, a sequence of matrix orthogonal polynomials satisfy a 3-term 
recurrence relation: 

tVnit) = An+lVn+l{t) + Bn{t)Vn{t) + A+P„_i(t), n > (1.19) 

where V-i{t) = 0, An are non-singular matrices and -B„ are Hermitian = En- 
Strikingly, as reviewed in [8] , for some families of MOPs, the Vn have been shown to be 
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eigenvectors of right second-order differential operators, that is to satisfy and equation 
of the form 

-^K(t)A2(t) + ^Vnit)A,{t) + Vn{t)Ao{t) = r„P„(t) (1.20) 

dt^ at 

where the differential coefficients A2, Ai and Aq are matrix polynomials independent 
of n, of degree not higher than 2, 1 and respectively and r„ are Hermitian matrices. 
When the MOPs Vn are functions of a discrete variable, as will be the case in our 
example, a discrete sum replaces the integral in the orthogonality relations f ll.lSp and 
one would wish to identify families of matrix polynomials that are eigenvectors of right- 
difference operators. The reader is referred to [H] for some specifics of matrix orthogonal 
polynomials of a discrete variable. 

1.5. Outline 

The paper is organized as follows. The recurrence relation that the matrix elements ilJn,k 
satisfy is obtained in Section [2l Upon analysis, it will be seen to imply that the ilJn,k are 
given in terms of 2 x 2 matrix orthogonal polynomials in k, Vn{k). The weight matrix 
W{k) for these MOPs is derived in Sect ion [3] from the unitarity of the representation and 
expressed using ipo^k and ipi,k solely. Special attention is given to those particular matrix 
elements in Section A. On the one hand, Hermite polynomials will be seen to enter in 
their explicit formula and on the other hand, they will be shown to form a /c-indexed 
doublet for which raising and lowering matrices will be constructed. This will provide the 
necessary tools to describe the bispectrality of the Vn{k) in Section [5] and to obtain the 
difference equation that they obey. Section [6] will be dedicated to the ladder relations 
and the Rodrigues' formula. It will be further shown in Section [7] that the matrix 
elements V'n.fc can be written as a "discrete" convolution involving Charlier and Meixner 
polynomials. Generating functions will be the object of Section [HI Introducing the 2- 
variable generating function of Perelomov and Popov [12] will bring to light the close 
connection between the matrix elements 'il)n,k and Hermite polynomials in two variables 
and will provide as a corollary intricate relations between the orthogonal polynomials 
involved. Finally the transformation law of the standard Hermite polynomials under an 
affine change of the independent variable is obtained as a by-product in Section [9] upon 
passing to the position representation. We close with concluding remarks and gather in 
an appendix the properties of orthogonal polynomials that are used in the text. 

2. Recurrence relations 

We shall initiate our analysis of the representation of Schi in the oscillator basis by 
determining the recurrence relation that the matrix elements ipn,k of (11.171) satisfy. We 
shall follow an approach already used in [T3] . 
We obviously have 

{k\a+aS\n) = k{k\S\n) = ktl)n,k (2.1) 
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and also 

{k\a'^aS\n) = {k\SS~^a~^aS\n) . 
With S given by (11.161) . using the Baker-Campbell-Haussdorf formula 

e^F e-^ = F + [X, Y] + ^[X[X, Y]] + ^[X[X[X, Y]]] + . . . 

we can straightforwardly obtain the automorphism of a~^a (a Schi algebra element) 
under S. In the following we shall use the four real parameters p, 6, a, 6 related to v 
and w as follows; 



(2.2) 



(2.3) 



iS if) 

v = ae , w = pe . 

We hence find 

S^^aS = a chp — e~*^ shp — ae^*'^, 
S-^a^S = a+ ch p - a e^^ sh p - a e^''; 

from where it follows that 

S-^a+aS = a+ach2p+ ^ch2p + a2 - ^ - la^sh2pe'' 



(2.4) 

(2.5a) 
(2.56) 



- a 
2^ 



sh 2 p e" 



+acx{ e*(^-^) sh p - e*"^ ch p) + a+a{ e''^^-^^ sh p - e"'^ ch p) . 



(2.6) 



Using (12. 6 p in (12. 2p and the actions of a and a+ on |n) given in (II. 9p . one finds 
upon combining (12.21) and (12.11) the following 5-term recurrence relation for the matrix 
elements ipn,k'- 



n + - ) ch2p + 



--\/n{n - 1) sh2pe*V«-2,/c ^ 



'Jn+2,k 



^ V(n + l)(n + 2) sh 2 p e' V. 
+v^(T(e*(^"^) shp- e^''chp)V'n~i,fc 

+V^rTTai e-*('-^) sh p - e"'^ ch p) (2.7) 

The analysis of this recurrence relation performed by setting n = 0, 1, 2, . . . and so on, 
shows that ip2n,k and ip2n+i,k are separately given in terms of ■?/'o,fc and ipi^k through 
expressions involving four sets of polynomials Pn{k), Qn{k), Pn{k) and Qn{k) of degree 
n in the variable k: 

V'2n,fc = Pn{k) tpO,k + Qn-li^l,k (2.8a) 
V'2n+l,fc = Pn{k) i'0,k + Qni'l,k- (2.86) 

Note that it is a polynomial of degree n — 1 that is the factor of ipi^k in (I2.8ap . The 
initial conditions are 



Po{k) =Qo(fc) = l 
g_i(A;) = Poik) = 0. 



(2.9a) 
(2.96) 



Observe also that in the two special cases p = and a = 0, the recurrence relation 
(12. 7p reduces to 3-term recurrence relations. In the latter case (p = 0), (12. 7p imphes in 
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particular that ipi^k is given in terms of V^o.fc- '^i,k = {cr — k/a) e^^ipi^k- These two special 
cases will be treated in Section [71 Until then and unless specified, we take p ^ a. 
When this is so, the two initial "values" ipo,k and ipi^k of the recurrence solutions are 
independent. 

In view of fl2.8al [27851) . it is natural to unify 4'2n,k and ip2n+i,k in a 2- vector that we 
shall denote "^nX- 

^n,k = ( "f'"'' I . (2.10) 
y 'ip2n+l,k J 

Let us now define the quantities ^n, Vn and Cn by 

e„ = -^y/{n-l)nsh2pe~'' (2.11a) 
T]n = V^a(e'(''-^)shp- e-'-^chp) (2.116) 
Cn= (n + ^)ch2p + a2-l (2.11c) 
and introduce the 2x2 matrices An and -B„: 

A„ = I ° 1 (2.12a) 

Bn = I 1 . (2.126) 

V2n+1 C2n+1 

Note that Bn is Hermitian B^ = Bn- 

It is then straightforward to see that the 5-term scalar recurrence relation (12. 7p can 
be written as the following 3-term recurrence relation for the 2- vector '^n,k'- 

k ^n,k = An+l^n+l,k + 5„ + (2.13) 

If we bring in the matrix polynomial 

/ Pn{k) Qn-l{k) \ 

Vn=[ ^ (2.14) 

it is also immediate to observe that the recurrence solution (I2.8al I2.86p -( l279al \2.9b\i now 
takes the following form 

^„,fe = Pn(A;)^o,fc (2.15) 

with initial condition 

Vo{k) = l2x2. (2.16) 

Evidently, f l2.15p and (12.161) could have been obtained directly as the recurrence 
solution of (I2.13p . Now remembering that for generic p and a, \E'o,A: has no annihilator 
(since otherwise ipo,k and ipi^k would not be independent), using (I2.15P in (I2.13P one 
arrives at the conclusion that the matrix polynomials Vn{k) obey themselves a 3-term 
recurrence relation, that is 

kVnik) = An+iVn+i{k) + 5„P„,fc(A;) + A+Vn-i{k), Vo{k) = 1. (2.17) 
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From the short introduction given in Subsection 11.41 one thus concludes that the 
polynomials Vn{k) are orthogonal. This is the central result of the paper, namely that 
the matrix elements f ll.l7p of the Schi representation are given in terms of the matrix 
orthogonal polynomials Vn{k) that are of degree n in the discrete variable k. Much of 
the following will be aimed at characterizing these Vn{k). 



3. Weight matrix 

We shall determine in this section the 2x2 weight matrix W{k) for which the polynomials 
Vn{k) defined by the recurrence relation fl2.17p are orthogonal. We shall see that it is 
readily provided by the unitarity of the representation. 
Since 5*+ = S~^, we have 

^n,m = {m\S~^S\n) (3.1) 

oo oo 

= Y,HS^\k) {k\S\n) = ^ iJn,ki^^,k- (3.2) 

fc=0 fc=0 

If we use the 2- vector '^n,k, this identity can be cast in the form 



^-.fc^ifc = K'nl2x2. (3.3) 

fc=0 

Using the expression fl2.15p giving '^n,k in terms of Vn{k), we see that (13. 3p implies 

oo 

J]Pn(/c)*0,fc^^fc^m(A;)+ = <5n,,nl2x2. (3.4) 
k=0 

This is the orthogonality relation we were looking for (see f ll.lSp in the continuous case), 
since it can be written 

oo 

^ Vn{k)W{k)V+{k) = Sn,ml2>c2, (3.5) 
fc=0 

where W{k) the 2x2 weight matrix is given by 



W{k) = ^o,kKk =1 I • (3.6) 

In the next section, we shall focus on the matrix elements ipQ^k and il'i,k, the initial values 
of the recurrence relations, that form as we see the elements of the weight matrix. 



4. The matrix elements 'j/'o,fe and "01,*; 

We wish to have explicit expressions for ipQ^k and il'i,k- To that end, we shall first 
compute the vacuum or ground state expectation values ipo^i = (0|S'|0). 
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In order to carry out the calculation of this quantity we need so-called "disentangling 
theorems", that is tools to factor the group element S as ordered products of 
exponentials solely of monomials in the creation and annihilation operators. For the 
first factor of S, this is easily achieved using the formula 

e^+^= e-[^'^l/2e^- e^ (4.1) 

known to be true when A and B commute with their commutator. One has thus 

The factorization of the second factor in S is more involved but fairly well known 



in view of its wide use. It is performed in details in ref. [T4j. In a nutshell, the approach 
therein is as follows. One lets 



Sl{X) 



(4.3a) 



and 



and postulate that 5*1 (A) = 5*2 (A). Differentiating both sides with respect to the 
parameter A, using the B-C-H formula (12. 3p and the posited identity, yields differential 
equations for the functions /(A), g{X) and h{X). Solving with the appropriate initial 
conditions gives the following identity when A = 1: 



= exp 

• exp 

• exp 



th p(a+)^ 



ln(ch p) ( a^a + 



e'^th pa^ , {w = pe''^). (4.4) 
We are now equipped to compute ipo,o- Using ( 14. 2p and f l4.4p . we readily see that 
V^o,o = {0\S\0) 



-.2/2/ 



(chp)-i/2^0|exp((Te*^a)exp 



1 



-id 



thp)(a+)2 |0) 



-.2/2 



(chp) 



'1/2 



exp(^-ia2e*(25-e)thp 



(4.5) 



4- 2- ^o,fc 

A recurrence relation for the matrix elements V'o.fc is easily obtained from observing that 
(A;|^a|0) = {k\SaS-^S\0) = 0. (4.6) 

In the usual fashion one finds that 

5a5"^ = chpa + e"*^shpa+ + a(e-^'^chp+ e^^''"^) shp), (4.7) 
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which leads to 

Vk + 1 ch p V'o.fc+i + e"*^ sh pVk ■ipo,k-i 

+a{ e-'^ ch p + e'^-^-^) sh p) ^o,k = (4.8) 

when employed in (14. 6p . This 3-term recurrence relation has a simple solution. Let 

1/1 \ '^/'^ 

V'o.fc = ( 2 ^''^ V ^'^ ^^'^^ 
with po = 1. Upon substituting in ( 14. Sp one finds 

Pk+i + 2kpk^i - 2spk = (4.10) 

with 

v2 

Comparing with the 3-term recurrence relation (lA.ip of the Hermite polynomials, we 
see that pk = H^iy) and hence that 

1/1 \ ^'/s 

^o,fc = ^(2e-^'thpj //fc(s)^o,o (4.12) 



where Hk{s) are Hermite polynomials of the variable s given in (14. lip and ip^^ is as in 
dMD. 

The matrix elements ifji^k can in turn be obtained from the relation 



(A;|a^|0) = VA^(A; + 1|^|0) = {k\SS-^aS\{}) . (4.13) 
Using (I2.5al) . it implies 

+ 1 V^o,fc+i = - e"^^ sh p - a e'^^V^cfc (4.14) 
in other words, that ipi^k can be expressed in terms of ipQ^k and ipo^k+V- 



V'l.fc = -- 7— [c^e-^%fc + l^o,fc+i]- (4.15) 
shp 

Recalling expression (I4.12p for ip^^k) one has for ipi^k'- 
1 1 fe-'^ , \^/2 



x[ae'^'-'^Hk{s) + (^e^nhp) i/fc+i(s)} ■ V'o.o- (4.16) 

Let us point out that this result could have been obtained using a different route. One 
could have started from 

V'l.fc = (fc|5a+|0) = {k\Sa+S~^S\Q). (4.17) 
Using the Hermitian conjugate of (14. 7p . this is seen to entail 
ifji^k = ch pVk i^o,k-i + e*^ shpVk+1 ipo,k+i 

+ {ae'^ dip + ae''^^-'^^ shp) V'o.fc- (4.18) 
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The 3-term recurrence relation fl4.8p for ipQ^k can then be called upon to eliminate V'o,fc-i 
in fl4.18p and recover (14.151) . 

The upshot of Subsections 14.21 and 14.31 is that ipo^k and ipi^k and hence the weight 
matrix W{k) (13.61) of the matrix polynomials Vn{k) are given in terms of Hermite 
polynomials in the variable s (a function of the group parameters). It is also manifest 
from (l4.12p and (I4.16P that ipo^k and ipi^k are independent when p 7^ 7^ a. 

4-4- Raising and lowering matrices 

For purposes that will be clear in the next two sections, it is relevant to have explicit 
ladder relations for ipo^k and ipi,k- As it turns out, the relations ensuing from the identities 



Vk + 1 ^i^k+i = {k\aS\t) = {k\SS-^aS\i), (4.19a) 
\^^^,k-l = {k\a+S\i) = {k\SS'^a+S\t), i = 0,l (4.196) 

allow to express 4'i^k±i in terms of ipo,k and ipi^k- Let us see how this goes. 

Equation (I4.14p already provides us with the desired result for ipo^k+i- Using (I2.5al) 
in (I4.19ap with i = 1, gives similarly: 

Vk + 1 V'l.fc+i = chpxfjo^k - (T(i~"^i'i,k - V2e"^^ shpV'a.fc. (4.20) 
Now, when n = 0, the 5-term recurrence relation (12. 7p gives ip2,k in terms of ■^o.fc and 

- a/2 ch p sh p e"*V2,fc = ~ ^ + ^ ~ ^0.*^ 

-cr(e*(^-^)shp- e"'^ chp)V^i,fc. (4.21) 

Of course, as noted in Section [21 if p = we see that ipi^k is given in terms of ifjQ^k and 
Tpo^k and ipi^k are thus dependent; similarly, if a = ip2,k is expressed only in terms of 
ipo^k and does not depend on ipi^k- Using (14.210 to substitute in (I4.20p for ?/'2,fc yields 
after simplification 

v/^TT^i,fc+i = -^{k + 1-^2) V^o,fc - (ye'^'-'^ thpV^i,fc. (4.22) 
ch p 

The same approach, using ( 12.5&P in ( ]4.19&P leads to 

\/k ipo^k-i = ch pipi^k - o- e'^ipQ^k (4.23a) 
for i = 0, and to 

ie 

Vktl^i^k-i = — f-(A;-a2)^o,fc-fTe<^-'^)cothp^i,fc (4.236) 
sh p 

for i = 1 (after using (I4.2ip ). 

At this point, it is useful to return to the formulation (I2.10p where the pair ifjQ^k and 
ijji^k is viewed as the 2-vector \E'o,fc- Introducing the 2x2 /c-dependent ladder matrices 

—(J ch p 



The Schrddinger group and matrix orthogonal polynomials 13 

it is checked that the relations f l4.23p and f l4.14p . (I4.22p respectively translate into the 
lowering and raising relations: 

v4fc*o,fc = v^^o,fc-i, A^o,fc = Vk + l ^o,fe+i- (4.25) 

It is also immediate to verify that 

Ak+iAk - Ak-iAk = 12x2- (4.26) 

5. Difference equation 

We derive in this section the finite-difference equation that the matrix orthogonal 
polynomials Vn{k) satisfy. Together with the 3-term recurrence relation (12.171) . this 
will express the bispectrality property of these Vn{k). Recall that fl2.17p followed from 
k'4'n,k = {k\a'^aS\n)\ expectedly, the difference equation will "dually" stem from 

nipn,k = {k\Sa^a\n) = {k\Sa~^aS''^S\n). (5.1) 

The transform Sa^aS^^ of a~^a under the automorphism generated by 5* is evaluated this 
time using (14. 7p and its hermitian conjugate. Given that result, following a procedure 
by now very familiar, we find from (15. ip : 



n 



'4'n,k = X{k) 1pn,k + Vk fi{k) 1pn,k-l + V A; + l/i(fc) i)n,k+l 



where 



+ ^k{k~l)u iJn,k-2 + \/lkTTj(kT^ViJn,k+2 (5.2) 

A(A;) = \{k) = ('^ + + ^) ch 2p + sh 2 p cos(25 - ^) - ^ (5.3a) 

1^ = a{ch2pe-'^ + sh2pe'^^-^^) (5.36) 

z/ =^e-^^sh2p. (5.3c) 



Equation (15.20 is easily transcribed in the matrix formalism where 



^n,k = I "f'"-' 1 = Vn{k)^0,k. (5.4) 



One obtains: 

X{k)Vn{k)^o,k + VkpVnik - l)^o,fc-i + Vk + lpVnik + l)^0,fc+l 

+ ^/k{k-l)uVn{k - 2)^o,k-2 

+ ^/ik + l)ik + 2)Wnik + 2)^o,fc+2 = l{n)Pn{k)-^o,k (5.5) 

where 
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Here one uses the lowering and raising matrices Ak and At of the last section and 
the ladder relations fl4.25p . to write (15 .Sp in a form where all the (matrix) terms in 
the equation act on the 2-vector \I/o,A: which can then be factored out since it has no 
annihilator when p 7^ 7^ a as already discussed. This then gives the following equation 
for the matrix orthogonal polynomials Vn{k): 

\{k)Vn{k) + ^iVn{k - l)Ak + JiVn{k + l)Ak 

+vVn{k - 2)Ak-iAk + yVn{k + 2)Ak+iAk = ^{n)Vn{k). (5.7) 

Introduce the shift operators T± acting on the left as follows on functions of / of the 
discrete variable k: 

f{k)T^ = f{k±l). (5.8) 

Let now be the right-handed matrix-difference operator defined by 

Ak = uT^Ak-iAk + uTlAk+iAk + fiT^Ak + fiT+Ak + A(A;), (5.9) 
we see that the matrix orthogonal polynomials Vn{k) satisfy the difference equation: 

Pnik)Ak=^in)Vnik). (5.10) 

•f— ^ ^ 

It is straightforward to check that the operator is hermitian A^ = A^. Recall that 

the inner product on the space of 2 x 2 matrix functions of k is defined by 

00 

{V{k), Q{k)) = Y,nk)W{k)Q+{k). (5.11) 

The hermitian conjugate of an operator X on that space is hence defined by {V, QX) = 

•f— -f— 

(VX~^, Q). Since X{k) is real, to prove that A^ = A^, it suffices to show that the 

two pairs of operators T_Ak and T+Ak as well as T^Ak-iAk and T^Ak+iAk are both 

hermitian conjugate pairs. Let us look at the first. We have 

00 

{V{k), Q{k)T_Ak) = Y,V{k)W{k){Q{k - 1)A)+ 

fe=0 

00 

= Y,v{k)mo,k'^tk^tQ{k-ir 

00 

= j2nk)Vk^o,k^t-M^-^)^ 

00 

= J2nk + l)VkTl^o,k+i'^tkQikV 



k=0 

00 



j2nk)T+Ak^o,kn,kQiky 



k=0 



= {V{x)T+Ak, Q{k)) (5.12) 
where we used (13. 6p . f l4.25p and (15. Sp . which confirms that 

(r_A.)+ = T+A. (5.13) 
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One demonstrates similarly that 

{T^Ak-iAkV = TlAk+iAk (5.14) 
to complete the proof that = A^. 

6. Ladder operators and Rodrigues' formula 

Right-handed lowering and raising operators Lk and Rk can further be obtained for the 
matrix polynomials Vn{k). They are found to act as follows: 

Vn{k)Lk = QnVn-i{k), (6.1a) 

Vn{k)Rk = Qn+lVn+l{k) (6.16) 

where 




e„ = , . (6.2) 



In order to determine the expressions of and R^, we respectively start from the 
identities 



{k\S\n - 2) = {k\Sa^\n) = {k\Sa^S-^S\n) (6.3a) 

and 

V(n + l)(n + 2) {k\S\n + 2) = {k\S{a+'f\n) = {k\S{a+fS-^S\n.) (6.36) 

We here limit ourselves to a brief summary of the derivation since the steps are 
analogous to those followed to obtain the difference equation. We appeal again to 
(14. 7p and its hermitian conjugate the compute Sa'^S~^ and S{a^)'^S^^. The results 
are then used in (16. 3p to obtain the corollary relations between 'ipn±2,k and ipn,k, 4'n,k±i 
and 'ipn,k±2- One then translates these relations in terms of the 2- vectors '^n,k- As in 
Section [5l one finally employs the ladder matrices Ak and Ak to get identities for the 
matrix polynomials Vn{k) by eliminating \E'o,fc- This brings one to find that : 

E = i(ch2p+ l)rMfe_iA + ^e2^^(ch2p - l)TlAk+iAk 

+a[ e*^(ch 2p + 1) + e-*(^-^) sh 2 p ]T_ A 
+a[e*(^+^)sh2p+ e-'^^~^^\c\i2p - l)]T+Ak 



+ 



e^^h2p(A; + cr^ + ^) +ych2p(e2'^+ e-2'(^-^); 



2/ 2 

+ y(e2*^- e-2*(^-^)) (6.4) 



and 



Lk = {RkV (6.5) 
where one should remember (15.13^ and (15.140 . 
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Note now from formula ( 16.161) . that the repeated action of Rk on the constant 
matrix l2x2 gives 

(^^e}jVn{k). (6.6) 
This is the Rodrigues' formula for the matrix polynomials Vn{k). 

7. Decomposition in terms of Charlier and Meixner polynomials 

In this section we shall exploit the product structure of S to show that classical 
scalar orthogonal polynomials in a discrete variables are building blocks of the matrix 
orthogonal polynomials Vn{k). The Charlier and Meixner polynomials will hence be 
seen to enter the picture. 

Define the matrix elements of the two factors in S (11.161) : 

Xm,fc = (fc|e("^V) (7.1) 
^n,m=(m|e[""^'-^('^^)'l/». (7.2) 

In terms of these we have 

^Pn,k = {k\S\n) = (A;|e("™^)e["'"'-^(°^)'l/» 

oo 

771=0 
OO 

= ^ Xm,k^n,m- (7.3) 

We shall now calculate the matrix elements Xm,k and (pn,m- Note that this will 
provide by the same token, the matrix elements of S when p = or cr = 0. Indeed it is 
obvious from (17.31) that 

i^n,k = Xn,k when w = 0, (p = 0) (7.4a) 

and that 

^n,k = ^n,k when V = 0, {a = 0). (7.46) 
These will again be determined from their recurrence relation. 

7.1. The matrix element Xn.k^ Charlier polynomials 

Using (I7.4ap . the recurrence relation for Xn,k is obtained from (12. 7p by letting p = 0. 
One observes that it then only has three terms and that it reads : 

kXn,k = in + a'^)Xn,k - Vn a e~'^Xn-i,k - Vn + 1 a e~'\ri+i,k- (7.5) 

An examination of (17. 5p shows naturally that its solutions are of the form Xn,k = 
Pn{x)xo,k where Pn{k) is a polynomial of degree n in the variable k. If we let 

P„(A;) = 4^(-ae-^^)-"P„(fc), (7.6) 
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we see that the normahzed polynomials Pn{,k) obey: 

Pn+i{k) = {k-n- a^)Pn{k) - na^Pn^^{k). (7.7) 

Comparing with flA.9p . we find that Pn{k) are normalized Charlier polynomials with 

a = (T^. Reverting to standard Charlier polynomials with the help of (lA.lOl) . we thus 
have 

Pnik) = {-a''rCn{k;a''), (7.8) 

and 

X^,k = ^iae''rC^ik;a')xo,k. (7.9) 



There remains to evaluate xo,k = {k\ e'^'^"''"^ |0) which is done directly with the help 
of the coherent states properties. From (14. ip . one has (similarly to (14. 2p ): 

Hence, 

gm-sa+|Q^ = e^^/^e^l -t;) = e-^^/^t;). (7.11) 
using successively (ll.llap and f ll.ll&p . One then obtains 

Xo,fc = (fc|e"1o)= e-^l^{k\-v) 

I ^' e-™/2 = ^^(-ae-*-')^ (7.12) 



^k\ Vk\ 

with the help of (ll.llcp and v = ae^^. Collecting gives the following formula for the 
matrix elements Xn,k'- 

\/n\k\ 

in terms of Charlier polynomials. 

7.2. The matrix elements tfn,k- Meixner polynomials 



As per ( 17.4&p . the recurrence relation for ipn^k is the special case of (12. 7p when a = 0, 
that is: 

k(Pn,k = (n + ^^ch2p-^ i^n,k - ^ \/n{n - 1) sh 2p e*Vn-2,fc 

-^^/{n + l){n + 2) sh2pe-'Vn+2,fc. (7.14) 

It has 3-terms with increments of 2. The elements with even and odd n are decoupled 
and the solutions are clearly of the form[|| 

y^2n,k = Pn{k)(po,k, (7.15a) 

^2n+l,k = Qn{k)(pi,k- (7.156) 

I In this subsection, as in the preceding ones, the use of P„ (and Qn) to denote polynomials is generic. 
The identifications do not carry from one section to another — we trust there should be no confusion. 
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Substitution in fl7.14p gives for P„ and Qn the following recursions: 



kPn{k) 



2n + - ) ch2p 



Pn{k) - -^/2n{2n-l)sh2pe''Pn-i{k) 



->y{2n + l){2n + 2) sh 2 p e-'^Pn+i{k) 



(7.16a) 



kQn{k) 



2n + - ]ch2p- - Qn{k) 



- ^ ^2n{2n + 1) sh 2 p e''Qn^i{k) 

- - v/(2n + 2)(2n + 3) sh 2 p e-'^^Qn+iik) . 



Let us first focus on Pn{k). We can normalize f l7.16al) by setting 

Pn{k) -- 



i=(-le--sh2p) P„(fc), 



to find that fl7.16al) reduces to 



Pn+l{k) 



k ( 1 

U + - 

2 V 4 

1 

~4 



ch2p 



4J 



Pn(A;) 



n n — 



sh2 2pP,_i(A;). 



With c, a parameter defined by 

,2 



so that 



c = th> 0<c<l 
1 + c 



ch2p 



sh^ 2p 



4c 



1-c' 

eq. f l7.18p can be rewritten as follows: 



:i-c)^ 



Pn^x{k) 



+ c(n + 1/2) 



- 1/2)- 



1-c) 



-P„_i(x) 



(7.166) 



(7.17) 



(7.18) 
(7.19) 
(7.20) 



(7.21) 



Comparing with ( ]A.15|) . we conclude that Pn{k) coincide with the normalized Meixner 
polynomials with x = k/2, (5 = 1/2 and c as is. Using (]A.16p . we thus have in terms of 
the standard Meixner polynomials: 



Pn{k) = 

Recalling that (|)„ = 

.2)^ ~ 22" n! 



c- 1 



,k 1 , 
M„(-;-,c). 



2 ^ . . . + n — 1), we verify that 



1 3 
2 

1 (2n)! 



(7.22) 



(7.23) 



Putting (17.171) . (I7.22p . and (I7.23P together, we get after some simplification (with the 
help of (17191) and (E20D), that 



Pn{k) 



1 



e*"^th>M„(-;-,th^p). 



k 1 



(7.24) 
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The polynomials of (]7.15&|) and f l7.16fe|) are obtained similarly. Their recurrence relation 
is normalized by taking 



Qnik) 



■e~'^ sh2p 



1 



V(2n + 1)! 



(7.25) 



it then becomes 
-k 



- n + - ch 2p + 



4 



Qn{k) 



—n[ n + 



1\ sh> 



Qn-i{k)- 



(7.26) 



Using the same parameter c as in f l7.19p and f l7.20p . after simple algebra, one casts fl7.26p 
in the form 

'k-l [n + (n + 3/2)c]- 



Qn+i{k) 



1 - c 



— n n + 



1 



2/1-0 



Qn{k) 



;Qn~i{k). 



(7.27) 



Comparing again with (lA.lSp . we realize that in this case the Qn{k) are normalized 
Meixner polynomials with x = {k — l)/2, /3 = 3/2 and c unchanged. Passing to the 
standard Meixner polynomials with the help of (lA.16p . we therefore have 

'3^ 



Qnik) 



where 



_ 3 5 

n 2 2 



- + n-l 
2 



1 (2n + l)! 

)2n 



Bringing the pieces together ( 07.251) . fl7.28p . 07.291) ) with c = th^p still, gives 



2" n! ' " V 2 2 

Keeping 07.15al I7.15&P in mind, there now only remains to compute ipo^k and (pi^k 



(7.28) 
(7.29) 

(7.30) 



fi,k = {k\ e 



[wa^-w{a+)^]/2\ 



^ = 0,1, (7.31) 

to complete the determination of the matrix element y:>n,k- Let us point out that ipo^k is 
the overlap of Perelomov SU(1, 1) coherent state [7] with k. We need to look at the two 
cases z = 0, 1 separately. Using the disentangling theorem (4.4) : 



[wa'^-w{a+f]/2\ 



|0) = {chp)~^/^exp\--e-''thp{a+f |0) 



(chp)-/2^ 



2 



n=0 



ni 



e-*'thp \2n). (7.32) 



Consequently, 



V5o,2fc = (chp) 

¥^0,2^+1 = 0. 



-1/2 Vm- 



k\ 



1 N ^ 



(7.33) 
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1 



--e-Mip(a+f |1) 



n=0 



n=0 

from where follows that 

<^l,2k = 0, 

V5i,2fc+i = (chp)"^/^ 



(7.34) 



(7.35) 



In summary, given fl7.15al I7.15&I) and in view of (17.331) and (I7.35p . the only non-zero 
ipn,k are the even-even and odd-odd elements. Their fomulas are found to be 



'^2n,2k 



trlVmM e-<"-.»|5g; . M,. (.; i. th= ,) (7.36.) 



2k 



k\n\ 



-1)'= V(2fc + l)!(2n-M)! (thp) 



k+n 



2"+^= A;!n! 
3 



(chp)3/2 



■M„(fc;-,thV) 
when combining (I7.24p with (I7.33P and with (I7.35p . 



(7.366) 



7.3. The convolution for the general matrix elements 

The decomposition for the general matrix elements is obtained by bringing together the 
results of the preceding subsections. From (17. 3p . (I7.33P and (I7.35P we have that 



n,k 



^2n,k 
^2n+l,fc 



E 



X2m,k ^2n,2m 
X2m+l,k ^2n+l,2m+l 



(7.37) 



Combining f l7.13p and (I7.36ap . we find after simple algebraic simplifications 

(_l)k ^k 



^n,k = Vn{k)'^0,k = 

oo 

■E 



m=0 



\fk\ 2"n! 

(-(t2 e^(25-e) thp)™ 
2'"m! 



^2/2 gj(ne-fc<5) ^^^n ^ 



\ 



V l^^^^''^2™+i(m,a2)M„(m;|;thV) ) 



(7.38) 
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7.4- '^o^k revisited 

It is instructive to examine how the expressions for ipo^k and il'i,k that we obtained 
in Section H] through recurrence relations can be recovered from the expansion ( I7.38P 
in terms of Charher and Meixner polynomials. As we shall see, it requires the use 
of formulas obtained by Gessel [I5], that are not widely known, for sums of Charlier 
polynomials. 

For notational ease, let 

t= e'^^^-'hhp. (7.39) 

Since Mq{x; f3,c) = 1, from (I7.38P we have 



0,k 



oo ^ 

■y— 

^ ml 



-ah 



1 1/2 



(chp) 



(7.40) 

m=0 ^ ~ ^ \ (ch p)3/2 '-'2m+l 1 

The well-known generating function (lA.lOp gives '^'^^(jCm{x; a)t'^ /ml. Here 
however, the degree of the polynomials in the sums are 2m or 2m + 1. Fortunately, 
the needed results have been worked out and are recorded in the Appendix. We shall 
have to treat the two components of f l7.40p separately. Consider first the top one. Using 
dAlll), we find 



-ae-'^)^e-" 



[k/2] 

. -ah/2 _ 

^ (1+ tf^-^ 



(chp) 1/2 

{-k)2, {-lY 



1=0 



il 



2a^ 



(7.41) 



where [k/2] denotes the integer part of k/2. 
Noting that 



-k)e = {-k){-k + 1) . . . (-A; + 2£ - 1) = (-1 



,21 



kl 



{k-2^)l 



and writing 



2a^ 



y/2i 



k-2t 



a 



tV2 



2^7v^ 



(7.42) 



(7.43) 



we see that fl7.4ip can be cast in the form 



ikS 



0,fc 



/H(chp)l/2 



-<72(l+t)/2 



k/2 



[k/2] 
1=0 



-11 



(.l{k-2l)l 



k-2t 



(7.44) 
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From the explicit expression ( IA.3P for the Hermite polynomials, we observe that the 
expression in square brackets in (17.441) is actually equal to Hk{s) where 

s = -^{t'/' + t-'/^) (7.45) 
v2 

is exactly the variable defined in (14. lip . Observing further that in terms of t, ?/^o,o given 
in (14.51) reads 



, .1 s fc/2 

e-'''{l) =(^e-''thp) , (7.47) 



and that 



2 J \2 

we see that (17.440 coincides with the expression for ipQ^k given in (I4.12p . 

Let us now look at the bottom component of (17.401) . Using flA.lSp we have in this 
case that 

V'l,^ = Dj, + Ek (7.48) 



where 

Uk = — ere e 



i5 



(chp)3/2 

Z^(l + t)2^-fc £! 1, 2a2y ^'■^^''^ 

1=0 ^ ' 



and 



?? (- ^o-i^U^-a^/a ere 



(chp)3/2 

correspond to the two terms occurring in the sum on the r.h.s. of f lA.lSp . Comparing 
(I7.49ap with fTTiT]) we see that 

Dk = ^i^o,k (7.50) 
chp 

and thus we readily have 
As for Ek, with 

(-*W2«-*) = (-l)-'(^3|3^ (7.52) 
and using this time the factorization 

f^)'-"-' . f ' (,53) 



2a^J \t^/^J \2a 



The Schrddinger group and matrix orthogonal polynomials 



23 



one arrives at 
1 



Eh 



_f^e-^'')'=e~"'(^+*^/2 I 1 



[(fc-l)/2] 



t \(fc-l)/2 



1=0 



(7.54) 



As before the term in brackets is identified as Hk-i{s) with the help of (IA.3I) . so that 



/fc! (chp)3/ 
Combining f l7.5ip and fl7.55p . we find 

^i,fc = Dk + Ek 

VfclchpV 2 



V2 



(7.55) 



^0,0- 



(7.56) 



In order to recover (14.161) . we use the recurrence relation ( lA.ll) to replace the Hk-i{s] 
by 



(7.57) 



In view of (17^ . 



-^t'/'Hk{s) + kHk-i{s) 
v2 



(7.58) 



from where it follows that 



shp 



-2j<5 



.e*(^-^)-|ai7fc(3)+(-) if,+i(s) 1^0,0 



(7.59) 



in accordance with fl4.16p . Let us record that from f l7.56p . using Hk-i{s), we have 
alternatively to f l4.16p : 



l.k 



fk\ chp 



thp\fc/2 



k e*"/^ 
th^/2 p 



(7.60) 
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8. Generating functions and 2-variable Hermite polynomials 

We shall discuss in this section, generating functions for the Schi group representation 
elements. We shall examine in turn a 2-variable and a 1-variable generating function. 
As we shall see below, in the former case, 2-variable Hermite polynomials will appear 
and hence be shown to have connections with the special functions that have arisen so 
far. 



8.1. 2-variable generating function 



Let X and y be two complex indeterminates. Following Perelomov and Popov [T2|, we 
define the 2-variable generating function G{x, y) as follows: 

r^f \ ^ x" y"" 
G{x,y) = > ^ ^n,k 

k,n=0 



X y 



k,n=0 



{k\S\n). 



Clearly, ^(x, y) can be looked at as the matrix element of S between the coherent states 
\x) and \y) : 

G{x,y) = {x\S\y). (8.2) 

We shall rely on the properties fll.lla[ ll.ll&t ll.llcp of these special states, to first 
evaluate G{x,y). Once again we shall use the disentangled form of 5* as provided by 
formulas (14. 2 p and (14. 4p . It is straightforward to see that 



{x\ e 



-vv/2 —vx 



{x + v\. 



We have thus 



{x\S\y) 



'Vv/2 p—vx 



■ exp 



X + v\ exp 



2 



ln(ch p) ( a~^a + 



■ exp 



U) 



As an intermediary step in the determination of the r.h.s. of (18. 4p , one needs 

1- 



[x + v\ exp 



— (In ch p) ( a^a 



\y) 



{x + v\i) {i\ exp — (Inch p) ^a^a + -j \m) {m\y) 



e,m=o 

OO 

E 

l,m=0 



1 



(X + v) 



exp — ( m H — ) ln(ch p) 



OO 

1 - th^ p)'/' 5^ ^[(1 - th^ pY^'ix + v)yr 



m=0 



[l-th^ pf'^ exp[(l - th^ pf/\x + v)yl 



15) 
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where we used f ll.llcp and employed the identity chp = (1 — th^ p) 

With the help of it is easy to complete the calculation to find that 



G(a;,i/) = (l-thV)'/'-exp(-i 



vv e 

2 



-ie 



exp 
1 



e'^^thp^^ + (1 -th^p)^/^x?/ 



thpt>' 



1/2, 



+- e"^ th py'^+ {-V - e"*^ th p w)x + (1 - th^ pY^^vy 



Now let 



.6) 



(8.7) 



zi =x; Z2 = y. 

Given that the exponent in the last term of (18.61) is quadratic in the indeterminates, we 
can rewrite G{x,y) in the form 

1 _ 1 
-vv 

2 2 



G{x, y) = il- th^ p)^/^ ■ exp ( - ^vv - ^ e''^ th p v"^ 



■ exp 



U,i=l,2 



ZiZj/2) 



where the 2x2 symmetric matrix [aij] has for its elements: 
ail = e~*^thp, = — e*^thp, 

ai2 = a2i = -(l-thV)'/', 



and where 



Cl 
C2 



-f(l + e^^thp), 

V 



ch p 



(8.9) 

(8.10a) 
(8.10&) 



Consulting formula (IA.18p . we conclude that G{x,y) is up to a factor, the generating 
function of the 2- variable Hermite polynomials if/£^„(ci, C2). We found indeed that 

00 _ I, „ 
G{x,y) = > r^ ^n^k 

fe,n=0 



(chp) 



-1/2 



■ exp 

°° -k „,n 



k,n=0 

from where one immediately reads that 

1 



1pn,k 



k,n 



(je'\l+ e'^thp) 



ere 



i5 



ch p 



f8.111 



(8.12) 



recalling (14. 5 p and putting v = ae*^. In other words, the Schi representation matrix 
elements ijjn,k are essentially 2-variable Hermite polynomials. As a consequence the 
relation 'i/n,k = 'Pn(fc)^o,A: (in view of (I4.12p and (I4.16P provides a relation between 2- 
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7(ifeT)T^'='2n+l(Cl,C2) 



1 N k/2 



( Hk{s) 



■Vn{k) 



1 

sh p 



\ 



1/2 



(8.13) 



with ci, C2 and s, respectively given by fl8.10ap -( l8.10fe|) and (14. lip . Analogously, (I7.38P 
entails an expansion of the 2-variable Hermite polynomials iffc,n(ci, C2) in terms of 
Charlier and Meixner polynomials. 



8.2. 1-variable generating function 

It is also instructive to consider the 1-variable generating function Fk[y) defined by 



Fk{y) = y^—^ipn,k = y^-^{k\S\n). 
^ Vn! V V ^' 



(8.14) 



It is obviously related as follows to the 2-variable generating function G{x,y): 



00 _i. 

X 



G{x,y) = J2^Fdy) 



(8.15) 



which can thus be looked at also, as the generating function for the Fk{y). Let us return 
to expression (18. 6p for G{x,y), this time, treating y as a. parameter and keeping only x 
as a variable. Take 

1 \ 1/2 

-e-'^'thp) X (8.16) 



and 



t vsh 2 p V2 



G(2,?/) = exp 



p^/^ a e'% + -e''th py- 



It is a simple matter to check that G{x, y) = G{z, y) can be written as follows: 

ie^^thm-2T 

■exp(2r2; - z^) ■ ipo^. 

With the help of (]A.2p . we recognize in (I8.18P the presence of the generating function 
exp(2rz — z"^) of the 1-variable Hermite polynomials. We therefore have 

iS 



ae 



G{x,y)=exp ^-^y + ^ e"^ th py"^ 
. ch p 2 



E 

vfc=0 



1/1 \ '^/^ 

-(-e-^hpj H,{r)x' 1 >o,o 



(8.19) 
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from where we find, given (18.151) . the following 1- variable generating function for the 
Schi representation matrix elements 



Fk{y) = E 



y" 



n=0 



n,k 



^V2 



-id 



thpj 



k/2 



■ exp 



ere 



y + -e''thpy^ ■ Hk{r) ■ ^Jofi- 
Lchp 2 



.20) 



Interestingly, it is expressed in terms of 1-variable Hermite polynomials in the 
variable r of (I8.17p . which can also be written 

Q-ie/2 

where s, as defined in fl4.1ip . is the argument of the Hermite polynomials in the 
expressions of V'o.A; and ijji^k (see Section Hj). We can now have a last additional look at 
these two fundamental matrix elements in relation with the generating function. 
Note that 



Evaluating f l8.20p at ?/ = 0, gives 



(8.22) 



1 



/MV2 



- e 



-ie . 



k/2 



thp) Hk{s) ipofi, 



precisely the expression (14.121) obtained in Section HI 
Similarly 



dy 







y 



/fc! V2 



1 \ k/2 

e'Hhp' 



ere 



-Hdr) + 



^ie/2 



-Hdr) 



_ ch p V2 sh p ch p dr J 

which is found to coincide with (I7.60p . using the Appell property (IA.4p of the Hermite 
polynomials. 



9. Position representation and afRne transformations of Hermite 
polynomials 

The matrix elements %pn,k of the Schi representation have been obtained quite generally 
from the structure of the group and in particular through an embedding of its generators 
in the Heisenberg (envelopping) algebra. When particular realizations are brought in, 
typically, more special functions become associated to the symmetry group and further 
relations with the special functions already present in the general matrix elements often 
result . We provide a brief illustration of that in this section by showing how formulas 
for the affine transformations of Hermite polynomials can be obtained in this fashion. 
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In the position representation, the creation and annihilation operators take the 
famihar form fll.6p and the vector states are reahzed by the wave functions ^n{x) = 
{x\n) given in fll.lOp . where Hermite polynomials interven^. Upon exponentiating 
the (generalized) vector fields arising in this picture, the Schrodinger group Schi is 
represented projectively as a transformation group on functions. Denote by $f (x) the 
transform of under 5* which we can write 

<l>^(x) = {x\S\n). (9.1) 

Parallely, it can be equated to 

oo 

{x\S\n) = ^(x|A;) {x\S\n) 

oo 
fc=0 

and therein lie identities for the Hermite polynomials that generalize somehow their 
generating function relation as well as that of the Laguerre polynomials (see below). 
We shall take v and w real, thereby focussing on affine transformations. When 

s = e = o, 

We already noted in ffLT2|) - ffLT3|) that a-a+ = \^d/dx and a'^-{a+f = 2xd/da; + l. 
In this case, S thus consists of dilations with parameter p, followed by translations with 
parameter \/2a. Standard Lie theory gives 

{x\S\n) = e^'^'^/'i^ ■ e^(^/^^+^/2)$^(3,) 

= ef/^<l>n[ef{x + V2(x)]. (9.4) 

Let us first consider separately the cases p = and a = 0. When p = 0, from (19. 2p . we 
have 

oo 

^n{x + V2a) = Y^ Xn,k^k{x) (9.5) 

k=0 

where Xn,k are the specializations of ilJn,k to p = given in fl7.13p . Using this formula in 
(19. Sp . together with the precise form (ll.lOp of yields 

oo / N ^ 

H4x + V2a) = iV2 ar e^^'^^+'^'Z^) ■ ^^^"(^^ ^')^fc(a:) (9.6) 

k=0 

which provides an expression for the transform of Hermite polynomials under 
translations ([IE], [E]). Note that (19. 6p can be viewed as a generalization of the 
generating function identity (IA.2P which is recovered from (19. 6p when n = since 

§ In this section |a;) is the eigenstate with real eigenvalue x of the position operator (a + )/V2 and 
should not be confused with a coherent state. 
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When cr = 0, we recall from Section [7] that ipn^k = V^n,fc and that only ip2n,2k and 
V^2n+i,2A;+i are non-zero. Eqs. (19 ■4p and (19.21) therefore amount to the pair of equations: 

oo 

e^/'^snle^x) = ^¥^2„,2fc$2fc(x) (9.7a) 

oo 

e^/'$2„+i( e^x) = J2 V2n+i,2k+i^2k+i{x). (9.76) 

fc=0 

Again, substituting for $„(a;) and using the explicit forms of the matrix elements ^ 
given m (ESD, we find from ( 1TO1 19:751) after some simplifications, the following two 
formulas 

if2„(e^x) = Mle-p/2e(^=''-i)^V2 th> 



oo / , x^, 



k=0 



i/2.+i(e-x) = (^!i±l)le-p/2e(^^^-)-^/^ii^ 



oo / , 



k=0 

which offer expressions for the transforms of even and odd Hermite polynomials under 
dilations of their argument. In this case, equations (I9.8a[ 19. can be looked at as 
generalizations of the generating functions relation for the Laguerre polynomials, keeping 
in mind the connection flA.SaMX.S&l) between the Hermite and Laguerre families. Indeed, 
when n = 0, it is straightforward to check that (9.8) reduces to flA.7p . 

Formulas for the transform of Hn{x) under the general affine coordinate 
transformation 

x' = eP{x + V2p) (9.9) 

are obtained by compounding (19. 6 p and (19. Sp . which is in keeping with the convolution 
(I7.38P found in Section [7] for the general matrix elements. While this is straightforward 
to achieve, the resulting expressions are a bit cumbersome and will not be recorded here. 
Let us rather point out that in terms of the matrix polynomials Vn{k), we have 



HoJx') 



1 



2"+i/2^(2n+l)! ' j k=0 

where x' is as in (19. 5p . This equation follows from the transcription of (19. 2p in vector 
notation 

<l>f„(x) \ ^ 

= E ^n,kM^) (9.11) 
^2n+l{^) / k=0 

and from using our fundamental result, namely that '^n,k = 'Pn{k)^k{x)- The reader 
will remember, as shown throughout the paper, that the MOPs Vn{k) can be obtained 
in a variety of ways: recurrence relation, difference equations, Rodrigues' formula etc. 
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10. Concluding remarks 

In summary, this paper has been concerned with the Schrodinger group Schi in one- 
dimension — a fundamental symmetry group of non-relativistic physics. Basically, we 
have achieved an explicit and thorough determination of the representations of that 
group on the natural Heisenberg algebra module. We have found that the matrix 
elements of this representation are expressed in terms of matrix orthogonal polynomials. 

Our study has provided in fact, what is probably the most simple example of matrix 
orthogonal polynomials in a physical context. In view of the efforts currently deployed 
to develop the theory of matrix orthogonal polynomials, we trust that this example will 
provide insightful illustration and offer a useful test bed. 

Of particular interest is also that we could present a full symmetry analysis; here 
again a rare case in dealing with MOPs. This has allowed to realize a systematic 
and exhaustive characterization (3-term recurrence relation, orthogonality relations, 
difference equation, ladder operators, Rodrigues' formula, generating functions and 
scalar contents). 

Let us mention in closing possible directions for further study. The 2-mode 
extension has an obvious physical interest. Examining the possible supersymmetric 
extension, q- deformation or parabosonic version would certainly be worthwhile. Finally, 
it would generally be desirable to explore other physical and algebraic contexts where 
matrix orthogonal polynomials will arise. 
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Appendix A. A compendium of formulas for orthogonal polynomials 

In order to make the paper self-contained and to specify our notations and conventions, 
we collect in this Appendix most formulas for orthogonal polynomials that are used in 
the body of the article. 

Appendix A.l. Hermite polynomials flS^ 

The Hermite polynomials Hn{x) can be defined by the following 3-terms recurrence 
relations: 



Hn+i{x) + 2nHn-i{x) - 2xHn{x) = 



(A.l) 
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with Hq{x) = 1. They have the generating function 

Hn{x)— = exp{2xz - z^), (A.2) 



n=0 



and the exphcit expression 



where [s] denotes the integer part of s. 
They enjoy the Appell property 

^Hn{x) = 2nHn-i{x). (A.4) 
da; 

The Hermite polynomials can be expressed in terms of Laguerre polynomials L'^{x) (see 
below) : 

H,^{x) = {-ir2'-n\L-'/\x'), (A.5a) 
H2n+i{x) = (-l)"22"+in!xLy2(a;2). (A.56) 

Appendix A.2. Laguerre polynomials fT^ 

The Laguerre polynomials L"{x), a > — 1 obey the 3-term recurrence relation 

{n + - {2n + a + 1 - x)L^(x) + (n + a)L'^_,{x) = (A.6) 

and have the generating function 

oo 

J2 L':{x)z- = (1 - z)—' exp (A.7) 



n=0 



Appendix A. 3. Charlier polynomials lT3f 

The Charlier polynomials c„(x;a), a > are polynomials of a discrete variable and 
satisfy the 3-term recurrence relation 

— xcn{x; a) = ac„+i(x; a) — {n + a)cn{x; a) + ncn-i{x; a). (A. 8) 

We shall also record the normalized recurrence relation obeyed by the monic version 
Pn{x) = + 0(n — 1) of the polynomials. Here 

xpn{x) = Pn+i{x) + (n + a)pn{x) + napn-i{x) (A. 9) 

where 

c„(x;a) = ) pn{x). (A.IO) 



A generating formula for the Charlier polynomials is 



J2cnix;a)-= e* 1-- . (A.ll) 



n=0 
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Using umbral calculus, Gessel [T5] has obtained a formula for a sum involving only the 
polynomials of even degree. In our conventions it reads: 

^ ^ ^n\\ 2 J ^ (l + tW~^ i\\ 2a) ^ ' 

n=0 £=0 ^ ' 

Gessel also proves a more general result that can be specialized to a sum involving only 
polynomials of odd degrees to give: 



n=0 

oo 



a;pn(a;) = Pn+i{x) H ^ H — — Pn-i[x) (A.15) 



Appendix A. 4- Meixner polynomials [Wf 

The Meixner polynomials Mn{x;(3,c), f3 > and < c < 1, are also polynomials of a 
discrete variable. Their 3-term recurrence relation is 

(c - l)xM„(x; /3; c) = c{n + f3)Mn+i{x- /3; c) 

-[n + (n + /3)c]M„(x; /3, c) + nM„_i(x; /3, c). (A.14) 

The normalized recurrence relation is 

n + in + (3)c , , n(n + f3 — l)c 
1 — c (1 ~ c) 

where 

1 /c — 1\ " 

Mn{x;/3,c) = -^(^^) Pn{x) (A.16) 
with = /?(/? + + 

Appendix A. 5. Multivariate Hermite polynomials fT8^ 

The Hermite polynomials in several variables were introduced by Appell and Kampe de 
Feriet in terms of generating functions. They are denoted by 

) (A.17) 

where Xi, i = 1, . . . ,m are the variables and rij the degrees of the polynomials in the 
variable Xj. 

With A = \ aij] a fixed positive definite symmetric square matrix, these polynomials 
are defined according to 



oo 



exp 



E a.- x.t, - -taA = Y: ±r^H4x) (A.18) 



2 V ^ ni\...nj.' 

, 4,j=l J ni=0 



when m = 1, we recover the standard Hermite polynomials in one variable for an = 2, 
We make use of the case m = 2 in Section [8] of the paper. Note finally that they satisfy 
a biorthogonality relation with a companion set of polynomials (see [IB])- 
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